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The self-similar renormalization group is used to obtain expressions for the spectrum of the 
Hamiltonian with the Yukawa potential. The critical screening parameter above which there are 
no bound states is also obtained by this method. The approach presented illustrates that one can 
achieve good accuracy without involving extensive numerical calculations, but invoking instead the 
renormalization-group techniques. 
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I. INTRODUCTION 

Renormalization-group techniques are widely used in quantum field theory, statistical mechanics, and solid-state 
physics. Their usage in atomic physics is less customary. The aim of the present paper is to show how renormalization- 
group ideas can be applied for calculating the spectra of quantum-mechanical Hamiltonians with realistic potentials. 

As a model for illustration we opt for a Hamiltonian with the Yukawa potential. This choice is based on the special 
role of this potential in different branches of physics. In plasma physics it is known as the Debye-Hiickel potential, in 
solid-state physics and atomic physics it is called the Thomas-Fermi or screened Coulomb potential, and in nuclear 
physics one terms it the Yukawa potential. Among recent important applications of this potential we may mention 
its use in the models describing metal-insulator transition [1]. 

The problem of finding the energy levels for the Yukawa potential has received a lot of attention in literature. Several 
approaches have been used for solving this problem: the simple variational procedure [2] , use of atomic orbitals with 
a set of fitting parameters [3] , Rayleigh-Schrodinger perturbation theory [4,5] , method of potential envelopes [6-9] , 
an iterative procedure [10], and different numerical calculations [11-14]. 

In this paper we demonstrate how the problem can be treated by employing the self-similar approximation theory 
[15-19] based on the ideas of the renormalization-group and dynamical theory. The outline of the paper is as follows. 
In Sec. II we sketch the main steps of the procedure using the self-similar renormalization group. We present only those 
formulas that are necessary for understanding the following calculations; all details and mathematical foundation can 
be found in Refs. [15-19]. In Sec. Ill we apply the approach to the Schrodinger equation with the Yukawa potential. 
In Sec. IV we obtain the sequence of renormalized energies. The convergence of this sequence is governed by control 
functions defined from the minimum of multiplier. In Sec.V the procedure is applied to calculating the renormalized 
critical screening parameter. Emphasize that we are presenting here an analytical method, not relying on heavy 
numerical calculations. Despite its analytical nature, the method gives quite good accuracy for the found critical 
screening parameter. 



II. SELF-SIMILAR RENORMALIZATION 



We give here the general sketch of the procedure [15-19], without specifying the nature of the functions involved. 
Suppose a function f{x) is defined by a complicated equation that is being solved approximately. Employing a 
perturbative algorithm, we may get a sequence of approximations, {-F/cj^o' sought function f{x). To make 

the sequence convergent, we incorporate into the approximations Fk ~ Fk{x, u^), with fc = 0, 1, 2, . . ., a set of control 
hmctions Uk = Uk{x), so that the sequence {fkix)}'^^^ of the approximations 

fk{x) = Fk{x,Uk{x)) (1) 
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be convergent. 

Make the change of variables by defining a function Xk{f) through the equation 

Fo{x,Uk{x)) = f, x = Xk{f), (2) 
in which / is the new variable. With this change of variables, Eq.(l) yields 

Vkif) = fkixkif)). (3) 

The transformation inverse to Eq.(3) is 

fk{x) = yk{Fo{x,Uk{x)). (4) 

Construct an approximation cascade {yk} by requiring the self-similarity relation 

yk+p{f) = VkiVpif))- (5) 

The trajectory {yk{f)}'kLo of this approximation cascade is, according to Eqs.(3) and (4), bijective to the sequence 
{fk{x)}kLo of approximations in Eq.(l). Embedding the approximation cascade into an approximation flow and 
integrating the evolution equation for the latter, we obtain the evolution integral 

fk+i fif 

in which fk = fk{x) is a given approximation, = fl{x) is a renormalized self-similar approximation, and 

d 

Vkif) = Pk+l{Xk,Uk) - Fk{Xk,Uk) + (Wfe+l - Uk)-g;^Fk{Xk,Uk) (7) 

is the cascade velocity, where Xk = Xk{f) and Uk = Uk{xk{f))- The right-hand side of Eq.(6), that is t^, is the 
minimal time necessary for reaching the renormalized approximation ,fk-f-i{x)- 

A fixed point Vkif) of the approximation cascade represents, by construction, the sought function f{x) which can 
be obtained from transformation (4). At the fixed point the cascade velocity Vk{f) ^ as fc ^ oo. If the cascade 
velocity is not zero exactly, but only approximately, then we have not an actual fixed point, but a quasifixed point. 
For instance, assuming that Vk{f) » and Fk+i « Fk, from Eq.(7) we have 

d 

{uk+i - Uk)^ — Fk{x, Uk) = 0. (8) 
dUk 

which is a quasifixed-point condition. 

The convergence of the approximation sequence {.fk{x)}^^Q is equivalent to the stability of the cascade trajectory 
{yk{f )}kLo- The stability of the latter can be analysed by defining the multipliers 

= ^jVkif) (9) 

and 

_ 5Fk{x,Ukix)) 

"^k^x) = — — (10) 

5Fo{x,Uk{x)) 

These multipliers are images of each other being related by the transformations 

Mfe(/) = mk{xk{f)), mk{x) = Hk{Fo{x,Uk{x))). (11) 
The trajectory is locally stable when 

\mk{x)\ < 1, \fxk{f)\ < 1. (12) 

The multipliers (9) and (10) describe the local stability, at the step k, of the cascade trajectory with respect to the 
variation of initial conditions. This type of local multipliers can be called quasilocal multipliers [19]. Another type of 
local multipliers defined as 

_ nikix) 
mk-i{x) 
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characterizes the local stability, at the step k, with respect to the variation of the point k — 1. The latter multipliers 

can be termed ultra-local. 

Recall that control functions arc introduced so that to provide the convergence of the approximation sequence 
{/fc(^)}feLo> tl^S't is, the stability of the cascade trajectory {yk{f )}kLo- This suggests a way for the practical definition 
of control functions. The local multiplier (10) may be written as 

mk{x) = Mk{x,Uk{x)) (13) 

with 

M,(.,n) = ^^M^/^^i^. (14) 

To produce the maximal stability of the cascade trajectory, i.e., the fastest convergence of the approximation sequence, 
for each fixed value of x, we have to require the minimum of the absolute value for the multiplier (13) with respect 
to the control function Uk{x). In other words, the principal of maximal stability is the condition for the minimum of 
the multiplier modulus, 

min|Mfe(x,u)| = \Mk{x,Uk{x))\. (15) 

u 

This condition gives us a constructive definition of control functions. 



III. YUKAWA POTENTIAL 



Now we illustrate how the scheme of Sec. II applies for calculating the eigenvalues of the radial Hamiltonian 

2m dr"^ 2mr'^ r 

with the Yukawa potential. 

It is convenient to pass to dimensionless quantities by scaling the above Hamiltonian and reducing it to the form 

^-—2d^ + ^2 -■ (16) 

Here r G [0,oo); I = 0,1,2,...; and a is a positive constant called the screening parameter. To return back to 
dimensional quantities one has to make the following substitutions: 

A 01 TT H 

r mAr, a — - , H 



Write Hamiltonian (16) as the sum H = Hq + AH, with the first term being the Hamiltonian 

1 d2 m + 1) u 

with a Coulomb-type potential, where u is yet unknown quantity which will latter generate control functions. 
Employing some variant of perturbation theory in powers of the perturbation 

AH = H-Ho = ^— , (18) 

r 

we may construct a sequence of approximate eigenvalues, Ek, and eigenfunctions, '^/s, respectively, 

£;fc^i?(f(a,u), T^.^vi'Hr-,^), (19) 

where fc = 0, 1, 2, . . . enumerates approximations, while n = 0, 1, 2, . . . and I = 0, 1, 2, . . . are the quantum numbers 
labelling the energy levels. For the initial approximation one has the spectrum 

V? 

^0 = -2{n + l + ir 
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and the wave function 
in which 



(n + 2Z + l)! 



1/2 1 / o... 



2ur \ I ur \ ^21+1 ' ^ur 



n+l+l\n+l+l) V n+l+l) " Vn+Z+iy' 



exp -— ^ i^M-f^ , (21) 



r(n + / + l)(-r) 

^nl.'; T(m. A- 1. A- ^\(r,.-m 



r(m + Z + l)(n - m)!m! 



TO=0 

is an associate Laguerre polynomial. 

To find the subsequent approximations, wc could use the Raylcigh Schrodingcr perturbation theory. However this 
would involve the following complication. The whole spectrum of the Hamiltonian (16) contains, in addition to discrete 
levels, the continuous part. Therefore, we would have to deal, besides the summation over discrete levels, with the 
integration over continuous ones. 

To avoid this complication, when we are interested only in discrete levels, we may employ the Dalgarno-Lewis 
perturbation theory [20]. Then, for the fc-approximation one writes 



p=i 

k 



V'fc =V'o + ^AV'p. (22) 



The first correction for the eigenvalues is 



AE, = {i;o,AHiJo), (23) 
and the first correction to the eigenfunction is a solution to the equation 

{Ho - Eo)Mi = i^Ei - AH)4,o. (24) 
Solving the Dalgarno equation (24), one may calculate the second correction to the eigenvalues, 

AE2 = {AiIjuAH^o), (25) 

and so on. 

The Dalgarno equation (24) is a nonhomogeneous equation whose solution can be written as the sum of the general 
solution to the corresponding homogeneous equation plus a particular solution to the nonhomogeneous equation. The 
solution to the homogeneous equation is, as is evident, proportional to ipo- So we may set 

A^i = CV'o + V, (26) 

with the proportionality constant C defined by the normalization condition {ipk,fpk) = 1, and the function ip being a 
particular solution to the nonhomogeneous equation 

(Ho - Eo)ip = {AEi - AH)iPo- (27) 

Prom the normalization condition {tpi,ipi) = 1, for the function = ^/iq + Atpi, omitting the second-order term, one 
has 

(V'o,A^i) = 0. (28) 

With Eq.(26), this gives 

C = -{^o,^). (29) 
Following the scheme described, with the notation 
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(3=^{n + l + l), (30) 



we obtain the first correction for the eigenvalues of bound states, 



where the integral 

n 



Ib 



\ + 1)\Jq I n \ i\ 



(n + 2/ • 

can be expressed as 

/3 (/J + 1)"+2;+2 

through the Jacobi polynomials 



1 " 

m=0 

having the properties 

(1) = c;:+„ (-1) = 

Another integral in Eq.(31) is 

7o = lim In = 1. 



(n — m)!rn! 



In this way, the first approximation Ex = Eq-\- AE-i becomes 

V? — 2ulf3 

2(n + Z + l)2' 



where 

1^ ~ (l + ^)2n+2i+2 '-^n+2J+l"-^n P 



For the ground-state level, when n = Z = and /3 = a/2u, Eq.(32) reduces to 



E, = -u-'[l-a\, (33) 



where the notation 

(2u + a) 



^-i-T^r:^ (34) 



is introduced. 

The ground state plays a special role defining, when it becomes zero, the critical screening parameter ac, above 
which there are no bound states. Therefore, in what follows we consider the ground state. 
Eq.(27), with 
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writes 



V'o- 



(35) 



(36) 



The solution to Eq.(36) must be a bounded function, \(p{r)\ < oo, for all r € [0, oo). 
Let us present ip as the product 



<p{r) = tpo{r)g{ur) 



in which the second factor satisfies the equation 



where p{r) is to be defined from Eq.(36), which yields 



dp _2„ 2r 

— = (1 — ar)2re exp 

dr u 



The latter equation gives 



p(r) 



ar^ + {a-l)(r + 



1 



-2r 



1 - cr 



1 



2u + a 



2u + a 



r exp I — 



2u + a 



with an integration constant Ci. The equation for g{r) results in 



g{ur) = aur + (1 — a) 



where C2 is an integration constant and 



2ur 



In(Mr) + Ei{—ar) 



2Ei{2ur) 



Jlur 



UT 



Ei{ar) 



-dx 



is the exponential-integral function. 

The boundness of ip requires that Ci = 0. The additive term containing a function proportional to ^0 should be 
omitted in ip, since such a term has already been included into Eq.(26). The latter means that we have to put C2 = 0. 
As a result we obtain 



<^(r) = ^/ue-'''' {2(j{urf + (1 - cr) [l - e""'' - 2ur\n{ur) + 2urEi{-ar)\ } . 
For the normalization constant in Eq.(29) we find 



C=--a + {l-a) 



a , 2a 

In — ; 7b 



2u + a 2u + a 



with the Eulcr constant = 0.577215. 

After finding function (26), we can calculate the second correction (25), which yields 



where 



AE2 = u^Ji- UJ2, 



Ji = lim J2{u), J2 = J2(/3), 

u—>0 



(37) 



(38) 



(39) 



with P = a/2u, and 

J^{u) = 2 {2Cr + 2ar2 + (1 - a) [l - e""''/" - 2r Inr + 2rEi (-^)] } e-^^'+''^'' dr. 
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Thus, for the second approximation for the energy we obtain 

E2=Ei+ U^Ji - UJ2, 

where 

a 1 — cj 



Ji 



2 1 + /3 



,^ ,cr(l + 3/3) ,^ ,2 



(l + 3/3 + /3^)/3 + 
(1 + /?)(! + 2/3) ^ 1 + 2/3 



IV. RENORMALIZED ENERGY 

From the results of the previous section we derive the following sequence of approximations for the energy: 

En = , 



Ei=Eo + u'^ 



{2u + ay 



E2 = El 



2u^ 4u^a 2u^{2u'^ + 5ua -2a + Sa^) 

~l~ 77^ '. To" ~l~ 



2 (2u + a)2 (2m + Q!)3 (u + a)(2u + a)3 

8u''(2u2 + 5Ma + 4a2) 16m'' , (2u + a)^ 

In ■ 



(u + a) {2u + a)^ {2u + a)^ 4u{u + a) ' 
For multiplier defined in Eq.(14) we have 



Mk{a,u) 



dEk{a,u) /dEo{a,u) 



Substituting here the derivatives 



dEi 

du 



du 

BEq _ 
du 

12u^ 



du 



= u — 



+ 



{2u + a)2 (2u + af ' 



following from the sequence (41), we find 

Mi{a,u) = 



8u^ - 4(2 - 3a)u2 - 6a{2 - a)u + 
(2u + a)3 



The control function u{a) = Ui{a) is to be defined from the principle of maximal stability (15). To this end, we, 
try the equation 

Mi{a,u) = 0, 

which gives 

Su^ - 4(2 - 3a)u^ - 6q(2 - a)u + = 0. 
This cubic equation has three roots of which we have to select a real one satisfying the asymptotic condition 



lim u{a) = 1. 

a— »0 



The latter implies that if the screening parameter tends to zero, so that the Yukawa potential transforms into the 
Coulomb one, then one must return to the exact solution known for the Coulomb potential. Really, under a ^ and 
u ^ 1, from the sequence in Eq.(41) it follows that iJfe — > — | for all k. Eq.(45), with this asymptotic condition, 
yields the control function 



/ X 1 " 2 / 3 



in which 

If = 

where 



3a V3(3 + 20a - 27a^) 

oj = arctan — — , 

^ 4 + 9a - 27a2 

and the upper value of the screening parameter, below which Eq.(45) possesses yet a solution, is 

10 + 7^ , , 

ao = — = 1.056306. (47) 

Thus, solution (46) exists only in the interval < a < ao. For a > ao, Eq.(45) has no real solutions satisfying the 
derived asymptotic condition. 

For a > ao we need to find a minimum of the multiplier (43), which is not necessarily zero. This can be done by 
solving the equation 

|^Mi(a,w)=0, (48) 

which results in the control function 

u{a) = -lja = 0.322876a. (49) 

One may notice that when Eq.(44) has a solution, then the procedure is similar to renormalizing perturbative 
terms by means of a variational optimization [21-26]. However, as is shown above, this equation not always possesses 
physically reasonable solutions. While the principle of maximal stability (15) always provides us with a solution 
defining a control function. Therefore, this principle is more general than the simple variational procedure. 

Substituting the found control function into the approximations in the sequence (41), we get the renormalized 
expressions 

efc(a) = Eh{a,Uk{a)). (50) 



For example, when a < ao, we have 



^^^"^ = " 2(2. + 3a) ' ^''^ 



2u(u-|-a)(2u + 3a)2 V2u + 3ay ^ Au{u + a) " ^ ' 

Respectively, for a > ao, we have to substitute the control function (49) into Eq.(41). 

The self-similar approximation for the energy is to be defined from the evolution integral (6). When no additional 
constraints are imposed, the minimal number of steps for reaching a quasifixed point is, clearly, one, = 1. In the 
interval < a < ao, the cascade velocity, given by Eq.(7), is 

^;i(/) = e2(a(/))-ei(a(/)), (53) 



e2(a) = -y 
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where the function a{f), according to Eq.(2), is defined by the equation 



Eoia, uia)) = -^u^ia) = / (54) 



resulting in a = q(/). For the evolution integral (6), we have 

df 
viif) 



(55) 



where ei = 61(0;) and et, = 63(0:). 

We calculated the values of the renormalized energies ei (a) , 62 (a) , and {a) as functions of the screening parameter 
a. To characterize the accuracy of these approximations, it is convenient to introduce the maximal percentage errors 



£fe = max 



efc(a) - e{a) 
\e{a)\ 



X 100%, 



where e{a) is an exact value of the energy. Notice that this definition of the maximal error has sense only when e{a) 
is not close to zero. In our case this definition works for < a < ao- For a > ao, when e{a) 0, it is possible to 
redefine the maximal error by shifting the definition of the energy by a constant [27]. 

The maximal percentage error, defined as is explained above, is 2% for 61(0;) and for 62(0;) and 62(0;) it is 1%. The 
multipliers (10), with a instead of x, satisfy the stability conditions of Sec. II. The stability of the procedure means 
its convergence which is also evident from the improvement of accuracy. 

V. CRITICAL SCREENING 

An important quantity characterizing the features of the Yukawa potential is the critical screening parameter, that 
is, such a value of the screening parameter a = ac above which there are no bound states. This critical parameter 
is defined by the condition e{ac) = 0. For each approximation e/j(a) for the ground-state energy there exists the 
corresponding critical parameter ak given by the equation 

efe(Q!fc) = 0. (56) 

From the approximations efe(a) obtained in the previous section, we find the sequence of approximations for the 
critical screening parameters: 

ai = 1, a2 = 1.0833. (57) 

To employ the self-similar renormalization for the sequence {ak}, as is described in Sec. II, we compose a sequence 
{afc(A)} of the partial sums 

k 

afc(A) = ^(ai-ai-i)Af-, (58) 

i=l 

in which fc > 1 and ao(A) = . As is clear from Eq.(58), 

ak = lini afe(A). (59) 

Then the sequence of approximations 

ai{X)=aiXP\ 

a2{X) = aiXP' + {a2 - ai)Af= 
can be renormalized in the way prescribed by Sec. II. First, we define the expansion function A(/) by the equation 

ai(A) =aiAPi =/, 
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which gives 



Writing the cascade velocity as 



«!(/) = <»2(A(/)) - <»,(A(/)) = (02 - <«,)A»"(/), 



we come to the evolution integral 



L,(X) (a2-ai)AP^(/) 
Notice that q;2(-^) Eq.(60) depends also on the parameters 



t*. 



so that we may write 



Integrating Eq.(60), we obtain 



p = pi>0, q=^-l>0, 
Pi 



"2(A) = o:2{\p,q)- 



1/9 



ai - q{a2 - ai)\m* 

To define p and g', consider the sequence {yk{f)} of the terms yk{f) = Oik{^{f))- Thus, we have 

yi{f) = f, 



(60) 



(61) 



(62) 



(63) 



y2(/) = / + (a2-ai)(£y^ . 

For the latter sequence we can find the multipliers defined as in Eq.(14). As is obvious, Mi = 1 and 

M2{X,p, g) = 1 + (a2 - ai)(l + q)XP'^. 



(64) 



The values of p and q are to be chosen so that to satisfy the principle of maximal stability (15), with p and q playing 
the role of control functions. Since, according to condition (59), we must put at the end A — » 1, we can consider the 
multiplier (64) for A ~ 1. In the case when A > 1, the minimum of IM2I from Eq.(64) is provided by g = 0. But if 
A < 0, then this minimum can occur at g = 1. Here we keep in mind that q, given by Eq.(61), is an integer and that 
the difference a2 — cti is positive in agreement with Eq.(57). In this way, from expression (63) we derive 



q;2(A,p, 0) = aiX^ exp 



a2 



ai 



if q = 0, and 



a2(A,p, 1) 



alXP 



ai — {a2 — ai)XPt* 



(65) 



(66) 



when q = I. 

The effective time t* has the meaning of the minimal number of steps providing the renormalization of , when 
A — > 1. If we put A ^ 1 in the evolution integral (60) before the integration, then for a| we would get ai + (a2 — oii)t*. 
From here we sec that t* = gives cti; one step, that is <* = 1, leads to 0:2; and t* = 2 results in 2a2 — ai- Therefore, 
the minimal number of steps necessary for getting a nontrivial renormalization is ^* = 2. 

Eqs. (65) and (66) show that, when A = 1, then a2 does not depend on p. Consequently, we may write 
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a2{l,q) = 0:2(1, J), g). 



(67) 



Putting in Eqs.(65) and (66) t* = 2 and A = 1, we obtain 



and, respectively. 



a2(l,0) = aiexp(2^^— ^ ) (68) 



«2(l,l)= o "'o • (69) 



3ai — 2a2 
As the final answer we set 

a; = i[a;(l,0) + a2(l,l)]. (70) 
Substituting into formulas (68) and (69) the numerical values from Eq.(57), we have 

a2(l,0) = 1.1813, 02(1, 1) = 1.1919. 

Thence, Eq.(70) yields 

a*2 = 1.1906. 

This value of the critical screening parameter coincides with the result of numerical integration [11]. 

In conclusion, we have applied the self-similar renormalization theory [15-19] to calculating the energy and the 
critical screening parameter for the Schrodinger equation with the Yukawa potential. The calculated values are in 
good agreement with the results of numerical computation. This demonstrates that renormalization-group techniques 
can be successfully employed for solving the Schrodinger equation, not only with simple anharmonic potentials [28,29], 
but also for more realistic cases, such as the Yukawa potential that is often met in different physical problems. To 
get an accurate value of the critical screening parameter, we have used a method analogous to the algebraic self- 
similar renormalization [30]. We paid here attention mainly to the ground-state level, although the procedure we 
demonstrated is applicable to excited levels as well, but calculations become a little more cumbersome. We hope 
that the results we obtained are sufficient for showing the usefulness of renormalization-group techniques in quantum 
mechanics. 
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